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With the increase of the lattice distortion, the orthorhombic manganites i?Mn03 (i?= La, 
Pr, Nd, Tb, and Ho) are known to undergo the phase transition from the layered A-type 
antiferromagnetic (AFM) state to the zigzag E-type AFM state. We consider the microscopic 
origin of this transition. Our approach consists of the two parts. First, we construct an 
effective lattice fermion model for the manganese Sd-bands and derive parameters of this 
model from the first-principles electronic structure calculations. Then, we solve this model 
in the Hartree-Fock approximation (HFA) and analyze the behavior of interatomic magnetic 
interactions. We argue that the nearest-neighbor interactions decrease with the distortion 
and at certain stage start to compete with the longer range (particularly, second- and third- 
neighbor) AFM interactions in the orthorhombic ab-plane, which lead to the formation 
of the E-phase. The origin of these interactions is closely related to the orbital ordering, 
which takes place in the distorted orthorhombic structure. The model is able to capture the 
experimental trend and explain why LaMnOa develops the A-type AFM order and why it 
tends to transform to the E-type AFM order in the more distorted compounds. Nevertheless, 
the quantitative agreement with the experimental data crucially depends on other factors, 
such as the magnetic polarization of the oxygen sites and the correlation interactions beyond 
HFA. 

KEYWORDS: undoped manganites, phase diagram, magnetic interactions, first-principles cal- 
culations, effective model 



1. Introduction 

For the long time LaMnOs was regarded as a prototypical example of parent (or undoped) 
manganites, where the strong Jahn- Teller distortion was believed to coexist with the (layered) 
A-type antiferromagnetic (AFM) state. "^"^^ The origin of this AFM state was one of the most 
disputed points about one decade ago, right after the new wave of interest to the phenomenon 
of the colossal magnetoresistance in the manganite compounds has just emerg ed.6-1^) Despite 
many differences in details, all theories of that period of time seemed to agree that the Jahn- 
Teller effect plays an important role in the alternating population of the Sx^— and 3y^— 
orbitals (Fig. 1), which is primary responsible for the directional anisotropy of interatomic 
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magnetic interactions underlying the A-type AFM phase. Indeed, simple considerations for 



Fig. 1. (Color online) (a): experimental phase diagram of i?Mn03 versus temperature and ionic 
radius of rare-earth elements (from ref. 14). Magnetic phases are denoted as paramagnetic (P), 
A-type AFM (A), spiral AFM (S), incommensurate (IC), and E-type AFM (E). (b) and (c): spin 
arrangement in the orthorhombic ab-plane, which takes place in the AFM phases of the A- and 
E-type, respectively, (d): alternating 3x^—r^ and 3y^—r^ orbitals and main magnetic interactions 
in the ab-plane, which are responsible for the relative stability of the A- and E-states. 

the superexchange (SE) interactions suggest that the alternating (antiferro) ordering of the 
3x^— and 3y^— orbitals in the orthorhombic ab-plane leads to the ferromagnetic (FM 
or F) coupling, while stacking (ferro) orbital ordering in the c-direction is responsible for the 
weak AFM coupling.^' ^^'^^^ 

The main surprise came later when it was found that after replacing La by smaller rare- 
earth elements (R), which systematically increases all kinds of the lattice distortions (including 
the Jahn- Teller one), the orthorhombic i^MnOa compounds undergo the change of the mag- 
netic ground state (Fig. 1).^^) Briefly, the least distorted LaMnOa forms the A-type AFM 
structure. The opposite-end compounds (starting from HoMnOs) form the so-called E-type 
(zigzag) AFM structure. In the intermediate region, the magnetic structure is incommensu- 
rate and keeps some features of the both A- and E-type AFM phases. The appearance of the 
E-type AFM structure, which breaks the inversion symmetry in otherwise centrosymmetric 
crystal environment, is particularly interesting. It can be hardly understood in terms of the 
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nearest-neighbor (NN) SE interactions alone, because such a mechanism would inevitably im- 
ply the change of the orbital state and operate against the large energy gain associated with 
the Jahn- Teller distortion. Therefore, it seems that the more realistic scenario should involve 
some longer range interactions. 

At the purely phenomenological level, the competition between the A- and E-type AFM 
phases in the ab-plane can be rationalized in terms of the following interaction parameters 
and trends (Fig. 1): 

• the NN interaction j|' , which, depending on its sign, favors either FM or bipartite AFM 
arrangement; 

• the 3rd-neighbor AFM interaction J3 which couples all 3rd-neighbor spins antiferromag- 
netically, as required for the E-type AFM structure. Therefore, J3 should be an indis- 
pensable ingredient of the model analysis. As we will see below, the main details of the 
magnetic phase diagram of RMnOs depend on the competition between j| and J3. If 
considered alone, the 3rd-neighbor AFM interactions would favor the formation of an 
infinitely degenerate group of states, including two zigzag AFM structures propagating 
along the orthorhombic a- and b-axes. The experimentally observed E-type AFM struc- 
ture is the one of them, which propagates along the a-axis and where the spins are 
antiferromagnetically coupled along the b-axis; 

• the 2nd- neighbor AFM interactions J^*, which lifts the degeneracy and together with J3 
determines the direction of propagation and the periodicity of the E-type AFM phase. 
The combination of and J3 appears to be sufficient to bind the directions of spins in 
each of the orbital sublattices, which are denoted as 3x'^—r^ or 3y'^—r^ in Fig. 1. 

Loosely speaking, if ferromagnetic j| dominates over J^" and J3, the magnetic ground state 
will be of the A-type. On the other hand, if the longer range interactions dominates, the 
magnetic ground state will tend to be of the E-type. The last ingredient, which stabilizes 
the E-type AFM phase is the small difference between the parameters j| acting in the FM 
and AFM bonds, which can be caused by either the exchange stiction or the orbital ordering 
effects. This difference is necessary in order to stabilize the directions of spins in two orbital 
sublattices relative to each other. 

The purpose of this work is to show that all these features are in fact closely related to the 
crystal distortion and the type of the orbital ordering realized in the orthorhombic i2Mn03 
compounds. We use the same strategy as in the previous work devoted to BiMn03.^^^ First, 
we derive an effective low-energy model for the Mn(3d) bands and extract parameters of this 
model from the first-principles electronic structure calculations based on the linear-muffin- 
tin-orbital (LMTO) method.^^) Then, we solve this model in the Hartree-Fock approximation 
and analyze behavior of interatomic magnetic interactions and the total energies. 
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Table I. Crystal structure parameters of i?Mn03 compounds, a, &, and c are the ortliorhombic lattice 
constants, Mn-0 are the interatomic distances, and ZMn-O-Mn are the bond angles (the first line 
is the angle in the c-direction and the second line is the angle in the ab-plane). All data are taken 
at room temperature except for LaMnOs, corresponding to 4.2 K. 





LaMnOa^^) 


PrMnOa^^) 


NdMnOs^o) 


TbMnOg^i) 


HoMnOs^^) 


a (A) 


5.532 


5.449 


5.416 


5.302 


5.257 


b 


5.742 


5.813 


5.849 


5.856 


5.835 


c 


7.668 


7.586 


7.543 


7.401 


7.361 


Mn-0 (A) 


1.906 


1.909 


1.905 


1.889 


1.905 




1.959 


1.953 


1.951 


1.946 


1.943 




2.188 


2.210 


2.227 


2.243 


2.222 


ZMn-O-Mn (°) 


157 


152 


150 


144 


142 




154 


151 


149 


146 


144 



The existence of the long-range magnetic interactions in LaMnOs was previously consid- 
ered in ref. 17, in the context of the local stability of the A-type AFM state with respect to 
other magnetic states. In the present work, we will further consolidate this idea and argue that 
it constitutes the basis for understanding the magnetic properties of all undoped manganites. 

The paper is organized as follows. In the next two sections we briefly discuss the main 
details of the experimental crystal structure (Sec. 2) and the electronic structure in the local- 
density approximation (LDA, Sec. 3). The construction of the model Hamiltonian for the 
Mn(3d) bands is considered in Sec. 4 and the strategy employed for the analysis of this 
Hamiltonian is briefly reviewed in Sec. 5. The behavior of interatomic magnetic interactions are 
discussed in Sec. 6. Sec. 7 is devoted to comparison with the experimental data. Particularly, 
we will consider the behavior of the correlation energies and the magnetic polarization of the 
oxygen sites, which is typically missing in the low-energy model. Finally, the brief summary 
will be given in Sec. 8. 

2. Crystal Structure 

All considered compounds crystallize in the highly distorted orthorhombic structure. The 
space group is -Dgfe Schonflies notations (No. 62 in International Tables). The primitive 
cell has four formula units. The crystal structure itself and its implications to the magnetic 
properties of LaMnOa have been discussed in many details in previous publications. Some 
crystal structure parameters are summarized in Table I. It also includes the references to 
the experimental lattice parameters, which have been used in the calculations. Generally, the 
crystal distortion in i^MnOs tends to increase in the direction La— ^^Pr^Nd^Tb^Ho. For 
example, such a tendency is clearly seen for the b/a and b/c ratios as well as for the Mn-0- 
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Mn angles. On the other hand, the Jahn-Teher distortion is not monotonous and takes the 
maximum in TbMnOs. For example, the ratio of the maximal and minimal Mn-0 bondlengths 
is 1.187 in TbMnOa (in comparison with 1.148 in the least distorted LaMnOa), and only 1.166 
in the following it HoMnOs. This structural anomaly is directly related to the anomaly of the 
crystal-field (CF) splitting, which will be discussed in Sec. 4. 

3. Electronic Structure in the Local-Density Approximation 

An example of the LDA band structure as obtained in the LMTO calculations for LaMnOa 
and HoMnOs is shown in Fig. 2. The LMTO bases, which was used in the valence part of the 



LaMnO 



HoMnO, 




-2 2 
Energy (eV) 



Fig. 2. (Color online) Total and partial densities of states as obtained in the local-density approx- 
imation for LaMnOa (left) and HoMnOa (right). The shaded area shows contributions of the 
manganese 3d states. Other symbols show the positions of the main bands. The Fermi level is at 
zero energy. 



spectrum, typically included the Mn(3d4sp), R{5d6sp), and 0{2sp) states. The i?(4/) states 
were treated as the (non-spin-polarized) core states. The atomic spheres radii were determined 
in two steps. First, we perform the LMTO calculations for the nominal composition, which 
includes 4 Mn, 4 R, and 12 O atoms, and find the atomic radii from the charge neutrality 
condition inside the spheres. Then, in order to better fill the unit cell volume and reduce 
the overlap between the atomic spheres, we add 12 to 16 empty spheres with the ls2p-basis. 
Typically, such a procedure guarantee a good agreement with the more accurate full-potential 
calculations. 

The electronic structure near the Fermi level is mainly formed by the Mn(3(i) states. There 
is also a considerable weight of the Mn(3d) states in the oxygen band. Due to the strong 
crystal-field (CF) effects in the MnOg octahedra, the electronic structure near the Fermi level 
splits into the "pseudocubic" Mn(eg) and Mn(t2g) bands. The Jahn- Teller distortion further 
splits the Mn{eg) band in two subbands lying at around 1 and 3 eV (Fig. 3). In NdMnOs, 
TbMnOs, and HoMnOs, these subbands are separated by an energy gap, whereas in the least 
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Fig. 3. (Color online) LDA energy bands for LaMnOs (left) and HoMnOs (right) as obtained in the 
original electronic structure calculations using the LMTO method and after the tight-binding (TB) 
parametrization using the downfolding method. Twelve low-lying bands spreading from around 
-1.0 till 0.4 eV are the "t2g bands" and the next eight bands are the "e^" bands. Notations of the 
high-symmetry points of the Brillouin zone are taken from ref. 26. 



distorted LaMnOs and PrMnOs, there is a small overlap between them. In the majority of the 
considered compounds, there is also a small overlap between upper Mn(eg) and R{5d) bands. 
An exception is HoMnOs, where these bands are separated by a small energy gap. 

4. Construction and Parameters of the Model Hamiltonian 

Our next goal is the construction of an effective model Hamiltonian for the Mn(3d) bands 
located near the Fermi level. For these purposes we use the method proposed in ref. 23. Many 
details can be found in the review article. ^^-^ The model itself is specified as follows: 

RR'«i«2 R {q} 

where c|^^ (cRa) creates (annihilates) an electron in the Wannier orbital centered at the 
Mn-site R, and a is a joint index, incorporating the spin (s= | or |) and orbital (m= xy, yz, 
z^, zx, or x^— y^) degrees of freedom. 

The one-electron Hamiltonian tRR' = ||tp|;^? || consists of the two parts: the site-diagonal 
elements (R=R') describe the crystal-field effects, whereas the off-diagonal elements (R^^R') 
stand for the transfer integrals, describing the kinetic energy of electrons. They are derived 
from the LDA band structure by using the formal downfolding method, which is totally equiv- 
alent to the use of the Wannier-basis in the projector-operator method. The comparison 
between the original LDA bands and the ones obtained in the downfolding method is shown 
in Fig. 3. In LaMnOs, the agreement is nearly perfect for the Mn(t2g) and the most of the 
Mn(eg) bands located in the low-energy part of the spectrum. In this region, the original elec- 
tronic structure of the LMTO method is well reproduced after the downfolding. Since upper 
Mn(eg) bands overlap with the La(5(i) bands, it is virtually impossible to reproduce all de- 
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tails of the electronic structure in the minimal model (1) limited to the five Wannicr-orbitals 
centered at each Mn-site. In this sense, the electronic structure obtained in the downfolding 
method is only an approximation to the original LDA band structure. Similar situation occurs 
in PrMnOa, NdMnOs, and TbMnOa. In HoMnOs, all Mn(3(i) bands are separated from the 
Ho(5(i) ones and well reproduced by the downfolding method. 

The one-electron parameters in the real space are obtained after the Fourier transforma- 
tion. Since we do not consider here the relativistic spin-orbit interaction, the matrix elements 
^RR'^ are diagonal with respect to the spin indices: i.e., t'^-^f =t^^^ Ssis2- Then, the site- 
diagonal part of tRR' = ||i]^R?^ II describes the CF effects. For example, the CF splitting is 
obtained after the diagonalization of ^rr. It is particularly strong for the Cg levels, being of 
the order of 1.5 eV (Fig. 4), and increases with the increase of the crystal distortion. As was 
pointed out in Sec. 2, some decrease of the e^-level splitting in HoMnOs in comparison with 
TbMnOs is related to the decrease of the Jahn- Teller distortion. For all considered compounds. 



•9 1 



ffi 



LaMnOj PrMnOj NdMnOj TbMnOj HoMnOj 



Fig. 4. Crystal- field splitting. Three low- lying levels are of the "t2g"-type and the next two levels are 
of the "eg"-type. 



the CF splitting is caused by the difference in the Mn(3d)-0(2p) hybridization in different 
Mn-0 bonds, which after the elimination of the 0(2p)-states gives rise to the site-diagonal 
elements in the model Hamiltonian. The effect of nonsphericity of the Madelung potential, 
which plays a crucial role in the t2g compounds,^^'^^) is relatively small for the eg-systems. 
For example in HoMnOs, it changes the e^-levels splitting by less than 3%. 

The directions of the CF splitting alternate on the perovskite lattice according to the 
D^f^ space group. The corresponding distribution of the eg-electron densities (or the orbital 
ordering) is shown in Fig. 5.^^^ As will be discussed in Sec. 5, this orbital ordering is directly 
responsible for the behavior of not only the NN but also the longer range magnetic interactions. 



7/26 



J. Phys. Soc. Jpn. Full Paper 




Fig. 5. Orbital ordering in LaMnOa derived from crystal- field eg orbitals of downfolded Hamiltonian 
(more specifically, the distribution of the electron density corresponding to the lowest Cg level in 
Fig. 4).^^^ Oxygen atoms are shown by small spheres. The vectors a, b, and c show the directions 
of orthorhombic axes. Other symbols show interatomic magnetic interactions in and between the 
planes, which are related to the given orbital ordering. 



Because of complexity of the transfer integrals, it is rather difficult to discuss the behav- 
ior of individual matrix elements of Nevertheless, some useful information can be 
obtained from the analysis of averaged parameters 

/ Xl/2 
f ,(^\ — I \ " fmim2fm.2mi 

\m1m2 / 

where d is the distance between the Mn-sites R and R'. All transfer integrals are well localized 
and practically restricted by the nearest neighbors at around 4A (Fig. 6). Already between 
the next nearest neighbors, the transfer integrals are considerably smaller. Generally, t^R' are 
larger for the least distorted LaMnOs and smaller for the more distorted HoMnOs. 

The screened Coulomb interactions 020304 have usual dependence on the spin indices: 
U^ia2a3a4=Um^m2m3m4^siS2^s3S4- Generally, the matrix C/^=||f/mim2m3m4 II can depend on the 
site-index R. The intersite matrix elements of U are considerably smaller. ■^^^ 

The matrix itself has been computed in two steps. ^^'■^^^ First, we perform the con- 
ventional constrained LDA (cLDA) calculations, and derive parameters of on-site Coulomb 
and exchange interactions between pseudoatomic Mn(3(i) orbitals. These parameters are typ- 
ically rather large because the do not include the so-called self-screening effects caused by the 
same 3d electrons, which participate in the formation of other bands due to the hybridiza- 
tion.24) The major contribution comes from the 0{2p) band, which has a large weight of the 
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Fig. 6. (Color online) Distance-dependence of averaged transfer integrals, iRR'(d) = 

.mim2j-m2mi"\ 1/2 
IZ^mima ''RR' ^R'R / 

Mn(3(i) states (Fig. 2). This channel of screening can be efficiently taken into account in the 
random-phase approximation (RPA) by starting from the interaction parameters obtained in 
cLDA and assuming that the latter already include all other channels of screening. ^^-^ All RPA 
calculations have been performed by starting from the LDA band structure. Nevertheless, in 
order to simulate the electronic structure close to the saturated (ferromagnetic) state, we 
used different Fermi levels for the majority (|-) and minority (J,-) spin states. Namely, it was 
assumed that the Mn(3(i) band is empty for the J,-spin channel and accommodates all 16 elec- 
trons (per one primitive unit) for the |-spin channel. Meanwhile, we get rid of the unphysical 
metallic screening by switching off all contributions to the RPA polarization function, which 
are associated with the transitions within the Mn(3(i) band.'^^^ 

Then, at each Mn site we obtain the 5x5x5x5 matrix of the screened Coulomb inter- 
actions. Since the RPA screening incorporates some effects of the local environment in solid, 
the symmetry of such matrices differs from the spherical one.^^-* Nevertheless, just for the ex- 
planatory purposes, we fit each matrix in terms of three parameters, which specify interactions 
between the 3d-electrons in the spherical environment: the Coulomb repulsion U=F^ , the in- 
traatomic exchange coupling J=(F^-|-F^)/14, and the "nonsphericity" i?=(9-F^— 5-F^)/441, 
where , F^, and F'^ are the radial Slater's integrals. These parameters have the following 
meaning: U is responsible for the charge stability of certain atomic configuration, while J 
and B are responsible for the first and second Hund rule, respectively. The results of such a 
fitting are shown in Table II. One can clearly see that the Coulomb repulsion U appears to be 
relatively small due to the self-screening effects, while J and B are much closer to the atomic 
limit. 

The model (1) does not explicitly include the oxygen states. This could be a serious prob- 
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Table II. Results of fitting of the effective Coulomb interactions in terms of three atomic parameters: 
the Coulomb repulsion U, the exchange coupling J and the nonsphericity B. All energies are 
measured in eV. 



compound 


U 


J 


B 


LaMnOs 


2.15 


0.85 


0.09 


PrMnOs 


2.07 


0.85 


0.09 


NdMnOs 


2.11 


0.85 


0.09 


TbMnOs 


2.24 


0.86 


0.09 


HoMnOs 


2.16 


0.85 


0.09 



lem in the case of manganites, which are known to be close to the charge-transfer regime. 
On the other hand, it is well know that in many cases a good semi-quantitative description 
of the magnetic properties of manganites can be achieved already in a minimal model com- 
prising only of the Mn(eg) bands. We will pursue the same point of view and concentrate 
on the behavior of the Mn(3d) bands. The magnetic polarization of the oxygen states will be 
considered in Sec. 7, where it will be also argued that this effect is partially compensated by 
correlation interactions in the Mn(3d) band beyond the Hartree-Fock approximation. 

5. Solution and Analysis of the Model 

The model Hamiltonian (1) was solved in the Hartree-Fock (HF) approximation."'^^' ^"^'^^^ 
After the solution for cacli magnetic state, the total energy changes corresponding to infinitesi- 
mal rotations of the spins magnetic moments near this state were mapped onto the Heisenberg 
model:^^'^^) 

RR' 

where or is the direction of the magnetic moment at the site R. The parameters {Jrr'} 
can be expressed through the one-electron (retarded) Green function, and the spin- 

dependent part of the onc-clectron potential, AVr, obtained from the self-consistent solution 
of the HF equations. For some applications, it is convenient to consider Jrr' as the function 
of the band filling: 

/U) 
duj'Jnn'i^'), (2) 
-oo 

where 

^RR'(^') = ^ImTVi. {gl^,{u;')AVn'gi,^{u;')AVR} (3) 

and Tri, is the trace over the orbital indices. In order to obtain the observable parameters, 
Jnn'i^) should be taken at the Fermi energy ep- Jur' = '^rr'(^f)- Some details of this 
procedure can be found in the review article^^^ as well as in the recent publication devoted to 
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6. Electronic Structure and Behavior of Interatomic Magnetic Interactions 

A typical example of the densities of states obtained in the HF calculations for the FM 
and several AFM phases of LaMnOa in shown in Fig. 7. Even in LaMnOa, which is the least 
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Fig. 7. (Color online) Densities of states obtained in the Hartree-Fock calculations for the ferromag- 
netic (F), A- and E-type antiferromagnetic phases of LaMnOa. The Fermi level is at zero energy 
(shown by dash-dotted line). Other symbols show the positions of the main bands. Different spin 
states are indicated by the arrows. 



distorted compound, the small value of [/, obtained in the combined cLDA-l-RPA approach, 
appears to be sufficient to open the gap in the eg band, so that all magnetic phases, including 
the FM one, become insulating. As expected, the increase of the number of the AFM bonds 
associated with the change of the magnetic state in the direction FM— >A-type AFM^E-type 
AFM results in the narrowing of all bands. Thus, the opening of the band gap is considerably 
facilitated by the interplay of the crystal distortion with the AFM arrangement of spins. For 
example, even small Jahn- Teller distorted appears to be sufficient to open the gap in the 
quasi-two-dimensional FM planes of the A-phase.'^^'^^^ A similar situation is expected for the 
quasi-one-dimensional spin chains in the case of the E-phase.^^^ In other compounds, with the 
increase of the crystal distortion the bandwidths will additionally decrease. In other respects, 
the position of the main bands is similar to the one displayed in Fig. 7. 

The distance-dependence of interatomic magnetic interactions Jrr' is shown in Fig. 8. 
One can clearly distinguish four types of interactions, which mainly contribute to the magnetic 
properties of i^MnOa: the NN interaction in the orthorhombic ab-plane, j|, which strongly 
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Fig. 8. (Color online) Distance-dependence of interatomic magnetic interactions, as obtained in the 
Hartree-Fock calculations for the ferromagnetic state. The interactions, which mainly contribute 
to the stability of the A- and E-type AFM phases, are shown in groups. The notations of these 
interactions are explained in Fig. 5. 



depends on the crystal distortion; the NN AFM interaction along the c-axis, J^; the 2nd- 
neighbor interaction in the ab-plane, J^, which operates along the orthorhombic b-axis; and 
the 3rd- neighbor AFM interaction in the ab-plane, J3, which operates only between those Mn- 
sites whose occupied Cg orbitals are pointed towards each other (see Fig. 5). Other interactions 
are considerably weaker. Particularly, the 2nd-neighbor interactions along the a-axis as well 
as the 3rd-neighbor interactions in the direction perpendicular to the occupied eg orbitals are 
small and can be neglected. 

The details of the behavior of the main magnetic interactions are shown in Fig. 9. The 
interaction j|' appears to be the most affected by the crystal distortion. When the crystal dis- 
tortion increases in the direction La^Pr^Nd^Tb^Ho, j|' gradually decreases and changes 
the sign at around Pr-Nd. Thus, the NN coupling in the ab-plane is FM at the beginning 
of the series and becomes AFM at the end of it. At the phenomenological level, such a be- 
havior can be related to the change of the orbital ordering in the Mn-O-Mn bond (Fig. 10). 
In LaMnOa, the Mn-O-Mn angle is closer to 180° (Table I). Therefore, the arrangement of 
the occupied Cg-orbitals at the neighboring Mn-sites is nearly "antiferromagnetic" which 
according to the Goodenough-Kanamori rules should correspond to the the FM coupling 
between the spins.^'^^'^^^ In HoMnOs, the deviation of the Mn-O-Mn angle from 180° is sub- 
stantially larger. Therefore, the "antiferromagnetic orbital ordering" is strongly distorted so 
that the spin coupling can become AFM. Nevertheless, as we will see below, although such a 
phenomenological interpretation is strongly affected by other details of the electronic struc- 
ture and particularly - by the hybridization between the t2g and eg states, which is caused by 
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Fig. 9. (Color online) The behavior of the main interatomic magnetic interactions for the i^MnOa 
compounds, as obtained in the Hartree-Fock calculations for the FM state: the nearest-neighbor 
interaction in the ab-plane, j| (a); the nearest-neighbor interaction between the planes, Ji (b); 
and the longer range interactions in the ab-plane, and J3 (correspondingly, b and c). The 
notations of the magnetic interactions are explained in figure 5. 




Fig. 10. (Color online) Fragment of the orbital ordering in the plane formed by the single Mn-O-Mn 
bond in the case of LaMnOs (left) and HoMnOa (right). 



the crystal distortion. 

Other magnetic interactions also depend on the crystal distortion. However, the distortion 
does not change the character of these interactions, and Jf*", J2 and J3 are AFM for all 
considered compounds. 

The most striking result of the present calculations is the existence of relatively strong 
longer range AFM interactions J2 and J3. The appearance of J3 is expected for the given type 
of the orbital ordering (Figs. 1 and 5). It operates between such 3rd neighbor sites R and R' 
in the ab-plane, whose occupied Cg orbitals are directed towards each other, and is mediated 
by the intermediate site, whose occupied eg orbital is nearly orthogonal to the bond (RR'). 
Although the direct transfer integrals between such sites R and R' are small (Fig. 6, note 
that the distance between 3rd neighbors in the ab-plane is about 8 A), the on-site Coulomb 
repulsion U is also relatively small (Table II). Therefore, the longer range AFM interactions, 
which are mediated by unoccupied Cg orbitals of intermediate Mn-sites, have the same origin 
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as the SE interactions, operating in the charge-transfer insulators via the oxygen states, ^^"^^^ 
and the mechanism itself can be called the "super-superexchange" . Another 3rd-neighbor 
interaction, operating between Mn-sites in the ab-plane whose occupied Cg orbital are nearly 
orthogonal to the bond connecting these sites, is negligibly small. A similar situation occurs in 
the low-temperature monoclinic phase of BiMnOs.^^^ The main difference is that the orbital 
ordering realized in BiMnOa is different from the one which takes place in the orthorhombic 
compounds. Therefore, the long-range AFM interactions in BiMnOa will tend to stabilize 
another magnetic state, which is also different from the E-state. 

The mechanism responsible for the appearance of the relatively strong interaction is 
not so straightforward. Nevertheless, as we will show below, some useful information can be 
gained from the analysis of the band-filling dependence of the 2nd-neighbor interactions in 
the ab-plane. 

Fig. 11 shows the behavior of the NN magnetic interactions as a function of the band 
filling. Somewhat unexpectedly, the NN interactions in LaMnOa are mainly formed by the 




Fig. 11. (Color online) Band- filling dependance of the nearest- neighbor magnetic interactions in the 
ab-plane (jj') and between the planes {Ji)- The magnetic interactions were calculated in the FM 
state for LaMnOa (left) and HoMnOa (right). Upper panel shows the behavior of the integrant (3), 
while the lower panel shows the exchange coupling (2). The Fermi level is at zero energy (shown 
by dash-dotted line). The positions of the t2g- and e^-bands are indicated by symbols. 
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t2g-band. Particularly, the values of both jj' and are well reproduced already after inte- 
gration over the t29-band spreading from -3.5 eV till -2.0 eV. The distribution of in 
the region of the occupied e^-band is antisymmetric. Therefore, there is a strong cancelation 
of contributions to coming from the bottom and the top of the occupied e^-band, so 

that the total integral (2) over the e^-band practically vanishes. In this sense, our explanation 
for the A-type AFM order in LaMnOs is rather different from the one adopted in the model 
calculations,^' "^^'"^"^^ which typically do not consider the rotations of the MnOg octahedra. Ac- 
cording to the present calculations, the behavior of the NN magnetic interactions in LaMnOs 
is mainly related to the hybridization between the atomic t2g- and e^-orbitals, which is induced 
by these rotations. Without the hybridization, all contribution of the half-filled t2y-band to 
the NN magnetic interactions are expected to be antiferromagnetic.^^' Our analysis shows 
that the hybridization can easily change the character of these interactions. 

The t2g-eg hybridization becomes even stronger in the more distorted HoMnOs, so that 
the contributions of the t23-band become ferromagnetic both for jj' and . On the contrary, 
all contributions of the e^-band to the NN interactions are antiferromagnetic. Therefore, the 
eg-band is totally responsible for the AFM character of NN magnetic interactions in the case 
HoMnOs. 

The behavior of 2nd-neighbor interactions in the ab-plane as a function of the band 
filling is shown in Fig. 12. Generally, the integrant ■yRR'(i^) oscillates in sign. Moreover, as 
the distance between the lattice centers R and R' increases, the number of such oscillations 
also increases. This property can be rigorously proven for the tight-binding bands, assuming 
that all transfer integrals (or "hoppings") are restricted by the nearest neighbors. Then, the 
number of nodes of ^rr'(i^) becomes proportional to the minimal number of hopes, which 
are required in order to reach the center R' starting from the center R.^*^'^^) Thus, ^rr'(i^) 
is expected to have more nodes for the 2nd-neighbor interactions in comparison with the NN 
ones, as it is clearly seen from the comparison of Figs. 11 and 12. Nevertheless, the lattice 
distortion and orbital ordering effects can cause some violation of these simple tight-binding 
rules. Let us consider the behavior of Jrr'(ci;) in the region of the e^-band, where jI{ijo) has 
only one node, which is qualitatively consistent with the tight-binding rules. Then, Jfiuj) 
has two nodes, which is again consistent with the tight-binding rules. Such a behavior is 
responsible for the strong cancelation of positive and negative contributions to Jf in the 
process of integration over uj and readily explains the fact that the final values of Jf are 
relatively small for all considered compounds. However, the cj-dependence of J2 appears to 
be strongly deformed. In the region of the e^-band it has only one node . Therefore, the strong 
cancelation, which took place for Jf , does not occurs for J^*. This leads to the strong anisotropy 
of the 2nd-neighbor interactions in the ab-plane, | J^'l '> | Jf I' which plays a vital role in the 
formation of the E-type AFM structure. Particularly, it readily explains the fact why the FM 
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Fig. 12. (Color online) Band-filling dependance of the second-neighbor magnetic interactions in the 
ab-plane. The magnetic interactions were calculated in the FM state for LaMnOs (left) and 
HoMnOa (right). Upper panel shows the behavior of the integrant (3), while the lower panel 
shows the exchange coupling (2). The notations of the magnetic interactions are explained in Fig. 
5. The Fermi level is at zero energy (shown by dash-dotted line). The positions of the t2g- and 
Cg-bands are indicated by symbols. 



zigzag chains in the observed E-type AFM structure propagate along the a-direction and are 
antiferromagnetically coupled along the b-axis (and not vice versa). 

Thus, the behavior of the main magnetic interactions replicates the gradual change of 
the crystal distortion. The form of both NN and long-range magnetic interactions is closely 
related to the orbital ordering realized in the distorted orthorhombic structure. Particularly, 
the crystal distortion explains 

• the gradual change of j\ from FM in the case of LaMnOa to AFM at the end of the series. 
Near the point of the FM-AFM crossover, j\ is small and the magnetic ground state is 
mainly controlled by the longer range interactions. 

• the existence of the longer range AFM interactions and J3, which bind the spin mag- 
netic moments within each orbital sublattice, and determine both the direction of prop- 
agation and the periodicity of the E-phase. 

Nevertheless, there should be an additional mechanism responsible for the relative orien- 
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tation of spin magnetic moments in two orbital sublattices, which are marked as 3x^— and 
3y2_^2 Yig_ 1. Since each spin in the E-type AFM structure participates in the formation 
of two FM and two AFM bonds with the nearest neighbors in the ab-plane, some difference 
between parameters j|' acting in the FM and AFM bonds is required in order to fix the di- 
rections of spins in the two orbitals sublattices relative to each other. Such a modulation 
of the parameters j|' can be caused by several mechanisms. Generally, once the symmetry is 
broken by the AFM spin order, orbital and lattice degrees of freedom will tend to adjust this 
symmetry change. 

One mechanism is purely electronic and related to the small deformation of the orbital 
ordering in the AFM phase. For example, in BiMnOa such a mechanism facilitates the for- 
mation of the tiit AFM structure, which breaks the inversion symmetry.^^) Nevertheless, in 
i?Mn03 the situation appears to be different. For all considered compounds, the NN inter- 
actions calculated in the E-phase satisfy the following condition: j|'(tt) < -/idi)) where the 
notations jj and jj, are referred to the FM and AFM bonds, respectively (Fig. 13). Thus, as 
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Fig. 13. (Color online) Nearest-neighbor magnetic interactions in the ab-plane of E-type antiferro- 
magnetic phase. The magnetic coupling in the FM and AFM bonds is denoted as tt and tJ,, 
respectively. 



far as the NN interactions are concerned, the E-phase appears to be unstable with respect 
to the spin rotations of two orbital sublattices relative to each other. Apparently, such a 
situation is realized in the intermediate region, corresponding to the IC- and S-states in Fig. 
1. Nevertheless, in order to stabilize the E-phase, we need another mechanism, which enforces 
the inequality j|'(tt) > </i(Ti)- Such a mechanism does exist and is related to the atomic 
displacements, which further minimize the total energy of the system via magneto-elastic 
interactions.^^' ^^-^ Although we do not consider it in the present work, from rather general 
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properties of the double exchange and SE interactions, it is reasonable to expect that the 
AFM character of j|(Ti) can be enforced by the conditions, which further enhance of the 
transfer integrals in the AFM bond.^''^ This can be achieved by either shrinking the Mn-Mn 
bond or increasing the Mn-O-Mn angle. The opposite distortions will favor the FM coupling, 
which are relevant to j|(TT)- 

7. Total Energies and Comparison with the Experimental Data 

In this section we consider the quantitative aspects of the problem. Particularly, we investi- 
gate whether the the experimental phase shown in Fig. 1 can be reproduced by the low-energy 
model (1) for the Mn(3(i) bands and, if not, which ingredients are missing in the model. 

We begin with the total energy calculations for the model (1) in the HF approximation 
(14). In LaMnOa, the lowest energy corresponds to the A-type AFM state, in agreement with 
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Fig. 14. (Color online) Total energies of different AFM states obtained for the model (1) in the 
Hartree-Fock approximation. All energies are measured relative to the FM states. The notations 
of the AFM states are standard for the manganites (see, for example, ref. 6). Apart from the A- 
and E-states, the C-state corresponds to the FM chains propagating along the c-axis, which are 
antiferromagnetically coupled in the ab-plane, and the G-state corresponds to the AFM coupling 
between all six nearest neighbors. 



the experiment. However, the next E-type AFM state is separated from the A-state by only 
1.1 meV per one formula unit. In PrMnOa and NdMnOa, the energy of the E-type AFM 
state appears to be lower than the one of the A-state, although experimentally both of these 
compounds are the A-type antiferromagnets (Fig. 1). Finally, for TbMnOs and HoMnOs, the 
model (1) yields the G-type AFM ground state, where all NN spins are coupled antiferro- 
magnetically. Thus, although the model (1) predicts the change of the magnetic ground state, 
it clearly overestimates the tendencies towards the antiferromagnetism, so that the transi- 
tion from the A- to E-type AFM state is expected in the wrong place (around PrMnOa and 
NdMnOs instead of HoMnOs). The correlation interactions beyond the HF approximation, 
will additionally stabilize the AFM states, and only worsen the agreement with the ex- 
perimental data. Therefore, before considering the correlation effects, one should find some 
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mechanism, which works in the opposite direction and additionally stabilizes the FM interac- 
tions. 

Such a mechanism can be related to the magnetic polarization of the oxygen sites. "^^'^^^ 
Although the model (1) is designed for the Mn(3(i) bands, the Wannier functions, which 
constitute the basis of the low-energy model (1), may have some tails spreading to the oxygen 
and other atomic sites. The weight of these tails in the Wannier functions is proportional to the 
weight of the 0(2p)-states in the total density of states for the Mn(3(i) bands (Fig. 2). In the 
case of the FM alignment of the Mn-spins, these tails will lead to some finite polarization at the 
intermediate oxygen sites (Fig. 15). Since the intraatomic exchange coupling Jq associated 




Mn O Mn 




Fig. 15. (Color online) Polarization of the oxygen sites caused by the tails of the Wannier functions 
centered at the manganese sites. In the perovskite structure, each oxygen site is located near the 
midpoint between two manganese sites. Then, in the case of the FM alignment, the tails from the 
Mn-sites have the same direction of spins, yielding the net magnetic moment also at the oxygen 
sites. In the case of the AFM arrangement, these tails cancel each other and the oxygen atoms 
remains nonmagnetic. 

with the oxygen atoms is exceptionally large, "^^"^^^ even small polarization can lead to the 
substantial energy gain. This contribution is missing in the model (1), where the form of 
the Coulomb and exchange interactions is assumed to be the same as in the limit of isolated 
Mn-atoms. In the case of the AFM alignment, the tails of the Wannier functions cancel each 
other and the net magnetic polarization at the oxygen sites is zero. 
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Table III. Magnetic polarization of the oxygen sites in different magnetic states of HoMnOs (namely, 
the absohitc values of the magnetic moments at the oxygen sites in fih). The first value was 
derived from the model analysis for the isolated Mn(3(i) bands, while the second value (shown 
in the parentheses) was obtained in the LSDA calculations, which also take into account the 
polarization of the 0(2p) band. Oab and Oc denote the oxygen sites located in the ab-plane and 
between the planes, respectively. Two lines in the case of the E-phase stand for the polarization 
in the FM (first line) and AFM (second line) Mn-O-Mn bonds. The finite polarization in some 
AFM Mn-O-Mn bonds is related to the oxygen displacements from the midpoint positions in the 
£>i? structure. 



phase 


Oab 


Oc 


F 


0.26 (0.11) 


0.23 (0.04) 


A 


0.25 (0.09) 


0(0) 


C 


0.08 (0.02) 


0.20 (0.03) 


G 


0.07 (0.03) 


0(0) 


E 


0.24 (0.07) 


0(0) 




0.09 (0.01) 



Below we present quantitative estimates of this effect for HoMnOs. By expanding the 
Wannier functions over the original LMTO basis functions,^^'^^^ one can find the distribution 
of the magnetic moments over all sites of the perovskite lattice in different magnetic structures. 
The obtained values of the magnetic moments at the oxygen sites, Mq, are given in Table 
III. The parameters Jq can be derived from the LMTO calculations in the local-spin-density 
approximation (LSDA).^^) It yields Jq = 2.1 and 2.2 eV for the oxygen sites located in the ab- 
plane and between the planes, respectively. Then, the energy gain, caused by the polarization 
of the oxygen sites, can be estimated from the formula AEq = —^Mq (with subsequent 
summation over all oxygen sites in the formula unit), which yields AEq = —102, —63, —29, 
—5, and —33 meV for the states F, A, C, G, and E, respectively. The effect is clearly too big. 
For example, by by combining these values with the total energies shown in Fig. 14, we would 
arrive to the FM ground state, which again contradicts to the experimental data. Then, what 
is missing? 

One effect is related to the polarization of the 0(2p)-band, which is not explicitly included 
in the model (1). It is true that since the 0(2p)-band is filled, it does not contribute to the 
total magnetic moment. However, it can contribute to the local moments, which cancel each 
other after the summation over the unit cell. Particularly, the polarization of the oxygen states 
in the 0(2p)-band appears to be the opposite to the one in the Mn(3d)-band, as it follows from 
the form of the Mn(3d)-0(2p) hybridization.^^) This effect is clearly seen by comparing the 
moments obtained for the isolated Mn(3d)-bands with results of the all-electron calculations. 



20/26 



J. Phys. Soc. Jpn. 



Full Paper 



Table IV. Correlation energies for several AFM states of HoMnOa measured in meV per one formula 
unit relative to the FM state. The correlation energies have been computed in the random-phase 
approximation (RPA), the second-order perturbation theory, and the T-matrix method starting 
from the Hartree-Fock approximation for each magnetic state. 



method 


A 


C 


G 


E 


RPA 


-4.9 


-19.5 


-24.7 


-14.1 


2nd order 


-6.7 


-14.9 


-17.9 


-10.3 


T-matrix 


-4.6 


-9.8 


-11.7 


-7.6 



which take into account the contributions of the 0(2p)-band (Table III). Indeed, the 0(2^?)- 
band substantiahy reduces the values of the magnetic moments associated with the oxygen 
sites (by factor two and more). Therefore, AEq will be also reduced. For example, by using the 
LSDA values for Mq (Table III), we find that AEq is reduced tih -13, -9, -1, -4, and -11 
meV per one formula unit for the states F, A, C, G, and E, respectively. By combining these 
AEq with the total energies shown in Fig. 14, we readily obtain the E-type AFM structure 
is realized as the ground state, in agreement with the experiment. The new values of the 
total energies, measured relative to the FM state, are —15, —15, —27, and —33 meV per one 
formula unit for the states A, C, G, and E, respectively. 

Another factor, which strongly affects the relative stability of different magnetic states, is 
the correlation interactions beyond the HF approximation. In order to estimate the energies 
of these correlation interactions, we tried three perturbative techniques starting from the HF 
solutions for each magnetic state. One is the random-phase approximation (RPA), which takes 
into account the lowest-order polarization processes, involving the excitation and subsequent 
deexcitation of an electron-hole pair.^^"^^) For these purposes, the RPA expression for the 
correlation energy has been adopted for the model calculations.^^^ Another method is the 
second order perturbation theory for the correlation interactions,^^' ^^'^^'^^^ and the third one 
is the T-matrix method,^^'^°) which takes into account higher-order effects. Results of these 
calculations for HoMnOa are shown in Table IV. Since the on-site Coulomb repulsion U is 
relatively small, all three methods provide rather consistent explanation for the behavior of 
the correlation energies, which tend to stabilize the AFM states relative to the FM one. The 
energy gain increases with the number of the AFM bonds in the direction F^A— >E— >C— >G. 
Thus, the correlation interactions act against the magnetic polarization of the oxygen sites and 
again tend to destabilize the E-state relative to the G-state. The situation is rather fragile and 
whether the E-state is realized as the ground state of HoMnOs depends on the delicate balance 
of these two effects and also on the approximations employed for the correlation energy. For 
example, RPA and the second-order perturbation theory seem to overestimate the correlation 
energy of the G-state and make the E-state unstable. On the other hand, the E-state, which 
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breaks the orthorhombic symmetry, should be additionaUy stabiUzed through the lattice 
relaxation. 

8. Summary and Conclusions 

On the basis of first-principles electronic structure calculations, we propose a microscopic 
model for the behavior of interatomic magnetic interactions in the series of orthorhombic 
manganites i^MnOs {R= La, Pr, Nd, Tb, and Ho), which explains the phase transition from 
the A-type AFM state to the E-state with the increase of the lattice distortion. Our picture 
is clearly different from the ones proposed in the previous studies. In fact, several authors 
emphasized the importance of the 2nd-neighbor interactions Jf and J2 in the orthorhombic 
ab-plane. For example, Kimura et al.^^^ considered the superexchange processes mediated by 
the 0{2p) orbitals in the distorted perovskite structure and argues that they can be responsible 
for the AFM interaction J2 and weakly FM interaction . Other authors^^' performed the 
mapping of the total energies derived from the first-principles electronic structure calculations 
onto the Heisenberg model and argued that under certain conditions Jf and J2 become 
comparable with j|' . However, such a mapping crucially depend on the form of the a priori 
postulated model, where the lack of some interactions (such as J3) can lead to an incomplete 
picture. In this sense, our approach to the problem is more consistent. 

• It does not make any a priori assumptions about the form of the Heisenberg model. 

• It goes beyond the conventional superexchange processes and takes into account other 
contributions to interatomic magnetic interactions.^^) 

Particularly, the contributions associated with the "supcr-supcrexchange" processes in the 
regime of relatively small on-site Coulomb interactions, give rise to the 3rd- neighbor coupling 
J3, which was overlooked in the previous studies. ^"^^ According to our point of view, J3 is one 
of the key players, which triggers the transition to the E-type AFM state in orthorhombic 
manganites. 

• The existence of J3 is directly related to the form of the orbital ordering. 

• J3 is responsible for the AFM coupling between 3rd-neighbor spins in the ab-plane, which 
is realized in the E-phase (Fig. 1). 

Since the longer range AFM interactions seem to be the intrinsic property of all undoped 
manganites, these interactions should be seen in the experiment, for example, on the inelastic 
neutron scattering. We expect the longer range interactions to take place even in LaMn03. 
Although it has A-type AFM ground state, the longer range interactions participate as the 
precursors of the E-phase, which is finally realized in the more distorted compounds. The 
neutron-scattering measurements on LaMn03 are available today. Nevertheless, the experi- 
mental data are typically interpreted only in terms of the NN interactions.^^' Definitely, 
the problem deserves further analysis. Particularly, in would be interesting to reinterpret the 
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experimental data by permitting the longer range interactions, particularly and J3. This 
point was already emphasized in ref. 17. It is possible that the longer range interactions are 
not particularly strong in LaMnOs, which has the highest Neel temperature (Tn, Fig. 1) 
and where the NN interactions clearly dominate. From this point of view, it would be more 
interesting to consider two other A- type AFM systems, PrMnOs and NdMnOs, which have 
smaller Tn and where the relative contribution of the longer range interactions to the magnon 
spectra is expected to be stronger. 

Although the proposed model is able to unveil the microscopic origin of the magnetic 
phase transition, the quantitative agreement with the experimental data crucially depends on 
the combination of the following three factors: 

• the correlation effects beyond the HP approximation; 

• the magnetic polarization of the oxygen sites; 

• the lattice relaxation in the E-phase, which breaks the inversion symmetry and gives rise 
to the multiferroic behavior. 

The detailed analysis of these effects presents and interesting and important problem for the 
future investigations. 
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